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1. Introduction 

The classical Fejer-Riesz factorization theorem gives the form of a nonnegative 
trigonometric polynomial on the real line, or, equivalently, a Laurent polynomial 
that is nonnegative on the unit circle. For the statement, we write D = {z : |z| < 1} 
and T = {^: |C| = 1} for the open unit disk and unit circle in the complex plane. 

Fejer-Riesz Theorem. A Laurent polynomial q{z) = X]feL_m 9'=^'^ which has com- 
plex coefficients and satisfies q(Q > for all £T can be written 

qiO = b(C)P, C e T, 

for some polynomial p{z) = Pq + piz + ■ ■ ■ + p,nz"^ , and p{z) can he chosen to have 
no zeros in D. 

The original sources are Fejer [22 and Riesz [17]. The proof is elementary 
and consists in showing that the roots of q[z) occur in pairs Zj and 1/zj with 
\zj\ > 1. Then the required polynomial p{z) is the product of the factors z — Zj 
adjusted by a suitable multiplicative constant c. Details appear in many places; 
see e.g. [H P- 20], [Ml P- 235], or [60l p. 26]. 

The Fejer-Riesz theorem arises naturally in spectral theory, the theory of 
orthogonal polynomials, prediction theory, moment problems, and systems and 
control theory. Applications often require generalizations to functions more general 
than Laurent polynomials, and, more than that, to functions whose values are 
matrices or operators on a Hilbert space. The spectral factorization problem is to 
write a given nonnegative matrix- or operator-valued function F on the unit circle 
in the form F = G*G where G has an analytic extension to the unit disk (in a 
suitably interpreted sense). The focal point of our survey is the special case of a 
Laurent polynomial with operator coefficients. 

The operator Fejer-Riesz theorem fTheorem l2.1[) obtains a conclusion similar 
to the classical result for Laurent polynomial whose coefficients are Hilbert space 
operators: if Qj, j — — m, . . . , m, are Hilbert space operators such that 

m 

Q(C) = E Q'^'^' C e T, (1.1) 

k— — m 

then there is a polynomial P{z) = Po + Piz + - ■ ■ + P„iz"^ with operator coefficients 
such that 

Q(C) = P(C)*P(C), C€T. (1.2) 
This was first proved in full generality in 1968 by Marvin Rosenblum [49; . The proof 
uses Toeplitz operators and a method of Lowdenslager, and it is a fine example 
of operator theory in the spirit of Paul Halmos. Rosenblum's proof is reproduced 
in gl 

Part of the fascination of the operator Fejer-Riesz theorem is that it can be 
stated in a purely algebraic way. The hypothesis (jl.l|) on Q{z) is equivalent to the 
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statement that an associated Toephtz matrix is nonnegative. The conclusion (jl.2p 
is equivalent to 2m + 1 nonlinear equations whose unknowns are the coefhcients 
Pq, Pi, . . . , Pm of P{z). Can it be that this system of equations can be solved by 
an algebraic procedure? The answer is, yes, and this is a recent development. The 
iterative procedure uses the notion of a Schur complement and is outlined in |j3l 

There is a surprising connection between Rosenblum's proof of the operator 
Fejer-Riesz theorem and spectral factorization. The problem of spectral factoriza- 
tion is formulated precisely in IjH using Hardy class notions. A scalar prototype 
is Szego's theorem (Theorem 14. 1|) on the representation of a positive integrable 
and log- integrable function w on the unit circle in the form \h\'^ for some func- 
tion h. The operator and matrix counterparts of Szego's theorem. Theorems 14.51 
and l4.7i have been known for many years and go back to fundamental work in the 
1940s and 1950s which was motivated by applications in prediction theory (see the 
historical notes at the end of Q. We present a proof that is new to the authors 
and we suspect not widely known. It is based on Theorem 14.31 which traces its 
origins to Rosenblum's implementation of the Lowdcnslager method. In i2]we also 
state without proof some special results that hold in the matrix case. 

The method of Schur complements points the way to an approach to multi- 
variable factorization problems, which is the subject of ^ Even in the scalar case, 
the obvious first ideas for multivariable generalizations of the Fejer-Riesz theorem 
are false by well-known examples. Part of the problem has to do with what one 
might think are natural restrictions on degrees. In fact, the restrictions on degrees 
are not so natural after all. When they are removed, we can prove a result, Theo- 
rem l5.H that can be viewed as a generalization of the operator Fejer-Riesz theorem 
in the strictly positive case. We also look at the problem of outer factorization, at 
least in some restricted settings. 

In recent years there has been increasing interest in noncommutative function 
theory, especially in the context of functions of freely noncommuting variables. 
In §6 we consider noncommutative analogues of the d-torus, and corresponding 
notions of nonnegative trigonometric polynomials. In the freely noncommutative 
setting, there is a very nice version of the Fejer-Riesz theorem (Theorem 16. 1|) . In a 
somewhat more general noncommutative setting, which also happens to cover the 
commutative case as well, we have a version of Theorem 15.11 for strictly positive 
polynomials (Theorem 16. 2[) . 

Our survey does not aim for completeness in any area. In particular, our 
bibliography represents only a selection from the literature. The authors regret 
and apologize for omissions. 
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2. The operator Fejer-Riesz theorem 



In this section we give the proof of the operator Fejer-Riesz theorem by Rosenblum 
|49j . The general theorem had precursors. A finite-dimensional version was given 
by Rosenblatt [48], an infinite-dimensional special case by Gohberg [26] , 

We follow standard conventions for Hilbert spaces and operators. If A is 
an operator, A* is its adjoint. Norms of vectors and operators are written || • ||. 
Except where noted, no assumption is made on the dimension of a Hilbert space, 
and nonseparable Hilbert spaces are allowed. 

Theorem 2.1 (Operator Fejer-Riesz Theorem). Let Q{z) — YTk^-mQ^z^ be a 
Laurent polynomial with coefficients in £(©) for some Hilbert space ©. If Q{C) > 
for all C G TT, then 

Q(C) = P(C)*P(C), CeT, (2.1) 

for some polynomial P{z) = Pq + Piz + • • • + Pmz"^ with coefficients in £(©). The 
polynomial P{z) can he chosen to he outer. 

The definition of an outer polynomial will be given later; in the scalar case, 
a polynomial is outer if and only if it has no zeros in D. 

The proof uses (unilateral) shift and Toeplitz operators (see [ITj and P^l)- 
By a shift operator here we mean an isometry on a Hilbert space such that 
the unitary component of S in its Wold decomposition is trivial. With natural 
identifications, we can write i3 = © © 6 © • • • for some Hilbert space © and 



S'(/io, hi, . 



(0, ho, hi, 



when the elements of are written in sequence form. Suppose that such a shift S 
is chosen and fixed. If T, A e £(^), we say that T is Toeplitz if S*TS = T, and 
that A is analytic if AS = SA. An analytic operator A is said to be outer if ran A 
is a subspace of i5 of the form ^ (B^ ® ■ ■ ■ for some closed subspace of 6. 

As block operator matrices, Toeplitz and analytic operators have the forms 



/To 


T-i 


T-2 






fAo 








...X 


Ti 


To 


T-i 






Ai 


Ao 













, A = 










Ti 


To 






A2 


Ai 






V : 










\ '■■ 






■■/ 



(2.2) 



Here 



' E*aS*'TEo\&, 
E^TS\^\Eo\e, 



J >0, 
J <0, 



(2.3) 



where Eog = ((7,0,0,...) is the natural embedding of © into fy. For examples, 
consider Laurent and analytic polynomials Q{z) — X^I-L-m Qfe-^'^ — 
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Pmz"^ with coefScients in £(©). Set Qj — for \j\ > m and 



Pj — for j > m. Then the formulas 



Tq 



/Qo 


Q-1 Q-2 


..A 




fPo 








..A 


Qi 


Qo Q-1 






Pi 


Po 













Tp = 










Q2 


Qi Qo 






P2 


Pi 


Po 




V 




••J 




\ '■■ 






■7 



(2.4) 



define bounded operators on f). Boundedness follows from the identity 

p CO 

/ (Q(C)/(C),5(C))ed^(C)= E (Q,-kh,gj)^, 



(2.5) 



where a is normahzed Lebesgue measure on T and /(C) = /o + /iC + /2C^ + ■ ■ • 
and g{() = go + giC + 52C^ + ■ ■ • have coefficients in ©, all but finitely many of 
which are zero. The operator Tn is Toeplitz, and Tp is analytic. Moreover, 



Q(C) > for all C e T if and only if Tq > 0; 

Q(C) = PiCyPiO for all C e TT if and only if Tq = T^Tp. 



Definition 2.2. We say that the polynomial P{z) is outer if the analytic Toeplitz 
operator A = Tp is outer. 

In view of the example (|2.4p . the main problem is to write a given nonnegative 
Toeplitz operator T in the form T — A* A, where A is analytic. We also want to 
know that if T = Tq for a Laurent polynomial Q, then we can choose A — Tp for 
an outer analytic polynomial of the same degree. Lemmas 12.31 and 12.41 reduce the 
problem to showing that a certain isometry is a shift operator. 

Lemma 2.3 (Lowdenslager's Criterion). Let f) he a Hilhert space, and let S G £(^)) 

be a shift operator. Let T G £(-$3) be Toeplitz relative to S as defined above, and 
suppose that T > 0. Let 9jt be the closure of the range ofT^^^ in the inner product 
of Sj. Then there is an isometry St mapping io^ into itself such that 

SrT'/^f = T^/'^Sf, fef). 

In order that T = A* A for some analytic operator A € £(^), it is necessary and 
sufficient that St is a shift operator. In this case, A can be chosen to be outer. 



Proof. The existence of the isometry 5*7- follows from the identity S*TS — T, 
which implies that T^^^Sf and T^^^f have the same norms for any f £ Sj. 

If St is a shift operator, we can view S)t as a direct sum Sjt ~ ®t ® ®t © ■ ■ • 
with ST{ho,hi, . . .) = (0, /iQ, /ii, • • ■ )■ Here dim(5T < dim6. To see this, notice 
that a short argument shows that T^I'^S^ and S*T^/^ agree on Sjt, and therefore 
T^/^(kerS'^) C kerS"*. The dimension inequality then follows because T^/^ is one- 
to-one on the closure of its range. Therefore we may choose an isometry V from 
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C5t into ©. Define an isometry W on Sjt into Sj by 

W{ho,hi,...)^ {Vh„,Vhi,...). 

Define A e £^{Sj) by mapping Sj into Sjt via T^/^ and then i^T into Sj via VF: 

A/ = / e jo. 

Straightforward arguments show that A is analytic, outer, and T — A* A. 

Conversely, suppose that T = A* A where A G S,{Sj) is analytic. Define an 
isometry W on S)t into by WT^/'^f = Af,fe S). Then WSt = SW, and hence 
S^"- = W*S*"'W for aU n > 1. Since the powers of S* tend strongly to zero, so do 
the powers of S^, and therefore St is a shift operator. □ 

Lemma 2.4. /n Lemma \2.3[ let T — Tq he given by (12. 4[) /or a Laurent polynomial 
Q{z) of degree m. If T = A* A where A S 'C(^) is analytic and outer, then A = Tp 
for some outer analytic polynomial P{z) of degree m. 

Proof Let Q{z) = Y.T=-m Qkz'"- Recall that Qj = for |j| > m. By dm apphed 
to A, what we must show is that S*-^ AEq — for all j > m. It is sufficient to show 
that S*"^~^'^AEo = 0. By applied to T, since T = A* A and A is analytic, 

E*A*S*^AEo ^ E^S*^TEo = Q.j =0, j> m. 

It follows that ranS'*™+M£;o L ra.nAS^'Eo for all /c > 0, and therefore 

ran S*"^^^ AEa -L ran A. (2.6) 

Since is outer, ran A reduces 5, and so ran S*"^'^^ AEq C ran A. Therefore 
S*'^+^AEo = by (HH), and the result follows. □ 

The proof of the operator Fejer-Riesz theorem is now easily completed. 

Proof of Theorem\2Jl Define T = Tq as in (gH). Lemmas O and [111 reduce the 
problem to showing that the isometry 5*7- is a shift operator. It is sufficient to show 
that ||S'^,"/I| ^ for every / in fiT- 

Claim: If / = T-^^^h where h ^ has the form h = {ho, . . . , hr,0, . . .), then 
S^^f = for all sufficiently large n. 

For if u G jj and n is any positive integer, then 

(^S^"f,T^/^u)^ =(^f,S'^T^/^u)^ ^ (T^/^h,T^/^S"u)^^ {Th,S"u)^. 

By the definition of T = Tq, Th has only a finite number of nonzero entries 
(depending on m and r), and the first n entries of 5"^ are zero (irrespective of u). 
The claim follows from the arbitrariness of u. 

In view of the claim, ||5'y"/|| for a dense set of vectors in S)t, and hence 
by approximation this holds for all / in 9)t- Thus St is a shift operator, and the 
result follows. □ 
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A more general result is proved in the original version of Theorem 12. II in [49j . 
There it is only required that Q{z)g is a Laurent polynomial for a dense set of g 
in (the degrees of these polynomials can be unbounded). We have omitted an 
accompanying uniqueness statement: the outer polynomial P{z) in Theorem 12.11 
can be chosen such that P(0) > 0, and then it is unique. See [2] and [50] . 



3. Method of Schur complements 

We outline now a completely different proof of the operator Fejer-Riesz theorem. 
The proof is due to Dritschel and Woerdeman [19^ and is based on the notion 
of a Schur complement. The procedure constructs the outer polynomial P{z) = 
Pq + Piz + • • • + PmZ™' one coefficient at a time. A somewhat different use of Schur 
complements in the operator Fejer-Riesz theorem appears in Dritschel TS]. The 
method in J8J plays a role in the multivariable theory, which is taken up in Sj^l 

We shall explain the main steps of the construction assuming the validity 
of two lemmas. Full details are given in [19] and also in the forthcoming book 
[3] by Bakonyi and Woerdeman. The authors thank Mihaly Bakonyi and Hugo 
Woerdeman for advance copies of key parts of [3] , which has been helpful for our 
exposition. The book 3J includes many additional results not discussed here. 

Definition 3.1. Let be a Hilbcrt space. Suppose T £ T > 0. Let ^ be a 

closed subspace of io, and let P^ £ be orthogonal projection of Sj onto M.. 

Then (see Appendix [XI Lemma there is a unique operator S £ £(.^), S > 0, 
such that 

(i) T--P*5P^^0; 

(ii) if S" e £(J^), S' > 0, and T - P^SPr > 0, then S <S. 

We write S — S{T, and call S the Schur complement of T supported on ^. 

Schur complements satisfy an inheritance property, namely, if C ^_|_ C Sj^ 
then S{T,n-) = S {S {T , Si+) , H^) . If T is specified in matrix form, 

T = ( ^ ^* ) : J? ® J?-^ -> J? © i^-" 



> 0. 



C 

then S = S{T,K) is the largest nonnegative operator in £,{R) such that 

'A-S B* 
B C 



The condition T > is equivalent to the existence of a contraction G £ £,{A,A-^) 
such that B = C^GA^ f Appendix \X[ Lemma [A. ip . In this case, G can be chosen 
so that it maps ran A into ranC and is zero on the orthogonal complement of 
ranj4, and then 

S = Ai{I-G*G)A^. 
When C is invertible, this reduces to the familiar formula S = A — B*C~^B. 
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Lemma 3.2. Let M e ii(io), M > 0, and suppose that 



M 



A B* 
B C 



for some closed suhspace ^ of f). 

(1) If S{M,A) = P*P and C = R*R for some P e £(il) and R e £(il^), then 
there is a unique X G such that 



M 



P* X' 
R* 



P 
X R 



and 



ranX C rani?. 



(3.1) 



(2) Conversely, if p.ip holds for some operators P,X,R, then S{M,K) = P*P. 
We omit the proof and refer the reader to [5] or [TOj for details. 

Proof of Theorem[K7\ using Schur complements. Let Q{z) — X]feL-m '5'='^'^ 
isfy (5(C) > for all C £ T. We shall recursively construct the coefficients of 
an outer polynomial P{z) = Po + Piz-\ h P^z" such that Q{C) = P{C)*P{C), 

C gt. 

Write ^ = © © © • ■ • and ©" © ■ • ■ © (5 with n summands. As before, 
set Qk — for |fc| > m, and define Tq e £^{S)) by 

/Qo Q-i Q-2 --A 



To = 



Qi Qo Q-1 
Q2 Qi Qo 

For each fc = 0, 1, 2, . . . , define 

5(fc) = 5(TQ,©'^-+i), 

which we interpret as the Schur complement of Tq on the first fc + 1 summands of 
i3 = 25 © • ■ • . Thus S{k) is a (fc + 1) x (fc+ 1) block operator matrix satisfying 

S{S{k),&+^) ^ S{j), 0<j<k<oo (3.2) 

by the inheritance property of Schur complements. 

Lemma 3.3. For each k = 0,1,2, ... , 

/ Yo {¥,■■■ Yk+i)\ 



5(fc + l) 

for some operators Yq,Yi, . 

(Yo Y, ■■ 



Y* 



S(k) 

V \Yk\J ) 
Yfc+i in £(0). For fc > TO - 1, 

Yk+i) = {Qo Q-i ■■■ Q-k-i) 
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Again see [3] or [T5] for details. Granting Lemmas |3 . 21 and we can proceed 
with the construction. 

Construction of Po,-Pi- Choose Pq = S'(O)^. Using Lemma [3T3l write 

In Lemma 15^ 1) take M — S{1) and use the factorizations 

5(5(1), ©1) ^ 5(0) = PJPo and 5(0) = P^Po . 
Choose Pi — X where X e £(Q3) is the operator produced by Lemma [3?2r i). Then 

^(l)^^? and ranPi CFiHPo. (3.3) 




Construction of P2 . Next use Lemma 

^'^> 

and apply Lemma [372^ 1) to M = 5(2) with the factorizations 

5(5(2),®^)^ 5(0) = Po*Po, 
5(1) = 

This yields operators Xi,X2 G ^{^) such that 



Po 


PA 








Po) 




Po) 



'P* X* X* 

5(2) =10 Po* Pi* I I Po 1, (3.4) 



(3.5) 




In fact, Xi = Pi. To see this, notice that we can rewrite (I3.4p as 

5(2) =P* 
^ ' yo R* ) \X Rj ' 

P=(t° n), X=(X2 Pi), R = Po 



Xi Po 

By p.3p and (|3.5p . ran Pi C ranPo and ranX2 C ranPo, and therefore ranX C 
ranPo. Hence by Lemma [3721 2). 

5(5(2), 02) = P*P=('^" ^Hf? (3-6) 
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Comparing this with 

S{S{2),&'') S{1) 





Pi) 




" ) 




Po) 




Po) 



(3.7) 



we get P*Pi = PqXi. By (IX5|), ranXi C ranPo, and therefore Xi = Pi. Now 
choose P2 = X2 to obtain 



f p* p* p* 
^0 ^1 ^2 



OOP 



ran 



, 



C ran 













Po 





\P2 


Pi 


Po 


Po 


') 




Pi 


Po 





(3.8) 
(3.9) 



Inductive step. We continue in the same way for all k = 1,2,3,.... At the fc-th 
stage, the procedure produces operators Pq, . . . , Pfc such that 



S{k) 



Pk\ (Po 



PoJ \Pk 







Po 



(3.10) 



ran 



C ran 



(3.11) 



By Lemma |3.3[ in the case k > m, 

( Qo {Q- 
/Q*-i\ 



S{k)^ 



Q- 



S(k-l) 



VV J 



(3.12) 



The zeros appear here when k > m, and their presence leads to the conclusion 
that Pfe = for A: > m. We set then 



P(z) -Po + Piz- 



Prr). ^ 



Comparing (|3.10p and (|3.12p in the case k = m, we deduce 2m + 1 relations which 
are equivalent to the identity 



Q(C) = P(C)*P(C), CeT. 
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Final step: P{z) is outer. Define Tp as in (|2.4p . With natural identifications, 

f Po {0 ■■■ )^ 

Tp = 




Tp 



(3.13) 



The relations (|3.1ip . combined with the fact that — for all k > ni, imply that 



ran 



2 C ran Tp 



V 



Hence for any G ©, a sequence /„ can be found such that 



Tpfn 



P2 



9- 



V 



Then by ((3T3)) . 



JnJ VO 

It follows that ran Tp contains every vector {Pog, 0,0, . . .) with g G £(35), and 
hence ran Tp 3 ran Pq (B ran Pq (B ■ ■ ■ ■ The reverse inclusion holds because by 
p.ll|) . the ranges of Pi, P2, ■ • ■ are all contained in ranPg. Thus P{z) is outer. □ 



4. Spectral factorization 

The problem of spectral factorization is to write a nonnegative operator-valued 
function F on the unit circle in the form F — G*G where G is analytic (in a sense 
made precise below). The terminology comes from prediction theory, where the 
nonnegative function F plays the role of a spectral density for a multidimensional 
stationary stochastic process. The problem may be viewed as a generalization of a 
classical theorem of Szego from Hardy class theory and the theory of orthogonal 
polynomials (see Hoffman [35l p. 56] and Szego [621 Chapter X]). 

We write and for the standard Hardy and Lebesgue spaces for the 
unit disk and unit circle. See Duren [5D]. Recall that a is normalized Lebesgue 
measure on the unit circle T. 

Theorem 4.1 (Szego's Theorem). Let w € satisfy w > a.e. on T and 

logw(C) da > — CX). 

Then w = a.e. on T for some h G H^, and h can be chosen to be an outer 
function. 
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Operator and matrix generalizations of Szego's theorem are stated in The- 
orems 14.51 and 14.71 below. Some vectorial function theory is needed to formulate 
these and other results. We assume familiarity with basic concepts but recall a few 
definitions. For details, see e.g. [30l[6T]) and [SOl Chapter 4]. 

In this section, 25 denotes a separable Hilbert space. Functions / and F on the 
unit circle with values in (5 and £(©), respectively, are called weakly measurable 
if {f{C),v) and {F{Qu,v) are measurable for all u,v G <3. Nontangential limits 
for analytic functions on the unit disk are taken in the strong (norm) topology for 
vector-valued functions, and in the strong operator topology for operator-valued 
functions. We fix notation as follows: 

(i) We write and ./^^(ig) for the standard Lebesgue spaces of weakly measur- 
able functions on the unit circle with values in © and 

(ii) Let and -ff^g) be the analogous Hardy classes of analytic functions on 
the unit disk. We identify elements of these spaces with their nontangen- 
tial boundary functions, and so the spaces may alternatively be viewed as 
subspaces of L| and ^^(g)- 

(iii) Let N^^^^ be the space of all analytic functions F on the unit disk such that 
(pF belongs to H^^^^ for some bounded scalar outer function ip. The elements 
of -^(i^g-) are also identified with their nontangential boundary functions. 

A function F e -^£(0) is called outer if FH^ is dense in for some closed 
subspace of ©. A function F G A^^jg) is outer if there is a bounded scalar outer 
function ip such that pF G -^^(g) and pF is outer in the sense just defined. The 
definition of an outer function given here is consistent with the previously defined 
notion for polynomials in SQl 

A function A G H^^^^ is called inner if multiplication by A on Sj^ is a partial 
isometry. In this case, the initial space of multiplication by A is a subspace of Sj^ 
of the form S)^ where 5^ is a closed subspace of (5. To prove this, notice that both 
the kernel of multiplication by A and the set on which it is isometric are invariant 
under multiplication by z. Therefore the initial space of multiplication by A is a 
reducing subspace for multiplication by z, and so it has the form Sj^ where 5^ is a 
closed subspace of © (see [29, p. 106] and [50l p. p. 96]). 

Every F G -^£(0) has an inner-outer factorization F = AG, where A is an 
inner function and G is an outer function. This factorization can be chosen such 
that the isometric set iJ| for multiplication by A on coincides with the closure 
of the range of multiplication by G. The inner-outer factorization is extended in 
an obvious way to functions F G N^^^y Details are given, for example, in [SUl 
Chapter 5]. 

The main problem of this section can now be interpreted more precisely: 
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Factorization Problem. Given a nonnegative weakly measurable function F on T, 
find a function G in -^£(0-) such that F — G*G a.e. on T. If such a function exists, 
we say that F is factorable. 

If a factorization exists, the factor G can be chosen to be outer by the inner- 
outer factorization. Moreover, an outer factor G can be chosen such that G(0) > 0, 
and then it is unique [50| p. 101]. By the definition of N'^^^y a necessary condition 
for F to be factorable is that 

/log+||F(C)|Ma<oo, (4.1) 
Jt 

where log''" x is zero or log x according asO<a;<lorl<a;<oo, and so we only 
need consider functions which satisfy (|4.ip . In fact, in proofs we can usually reduce 
to the bounded case by considering F/\(p\'^ for a suitable scalar outer function (p. 

The following result is another view of Lowdenslager's criterion, which we 
deduce from Lemma [2.31 A direct proof is given in |6T] pp. 201-203]. 

Lemma 4.2. Suppose F £ ^£(0) ci'i^d F > a.e. on T. Let be the closure of 
F's in Lg, and let Sp be the isometry multiplication by C on Kp. Then F is 
factorable if and only if Sp is a shift operator, that is, 

oo 

nC"^^={0}. (4.2) 

Proof. In Lemma 12.31 take = viewed as a subspace of , and let S be 
multiplication by ( on Sj. Define T £ ^{^) by Tf — PFf, f £ S), where P is 
the projection from onto H^. One sees easily that T is a nonnegative Toeplitz 
operator, and so we can define Sjp and an isometry St as in Lemma 12.31 In fact, 
St is unitarily equivalent to Sp via the natural isomorphism W : — > such 
that W{T^ f) = F^f for every / in f). Thus Sp is a shift operator if and only if 
St is a shift operator, and by Lemma [^31 this is the same as saying that T = A* A 
where A S £(i5) is analytic, or equivalently F is factorable [50l p. 110]. □ 

We obtain a very useful sufficient condition for factorability. 

Theorem 4.3. Suppose F G i^J®) o-nd F > a.e. For F to be factorable, it is 
sufficient that there exists a function tp in i^^^^j such that 

(i) i;F e if-^)/ 

(ii) for all ( £ T except at most a set of measure zero, ?/'((') |-F(^)© is one-to-one. 
If these conditions are met and F = G*G a.e. with G outer, then ipG* £ H'^f^igy 

Theorem 14.31 appears in Rosenblum [49] with '0(0 — C™ (viewed as an 
operator- valued function). The case of an arbitrary inner function was proved 
and applied in a variety of ways by Rosenblum and Rovnyak [501 [51] ■ V. I. Mat- 
saev first showed that more general functions -ip can be used. Matsaev's result is 
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evidently unpublished, but versions were given by D. Z. Arov 1, Lemma to The- 
orem 4] and A. S. Markus 41, Theorem 34.3 on p. 199]. Theorem 14.31 includes all 
of these versions. 

We do not know if the conditions (i) and (ii) in Theorem 14.31 are necessary 
for factorability. It is not hard to see that they are necessary in the simple cases 
dim © = 1 and dim © = 2 (for the latter case, one can use [SOj Example 1, p. 125]). 
The general case, however, is open. 

Proof of Theorem \4-3\ Let F satisfy (i) and (ii). Define a subspace Tl of by 

oo oo 

We show that 9Jl = {0}. By (i), 

OO OO OC 

.prim = i'Pi fl ("F^Hl C fl C'i^FHl C fl CHl = {0}. (4.3) 

n— n— n— 

Thus tpF^m = {0}. Now if 5 G an, then ^pF^g = a.e. by g3|). Hence pig^Q 

a.e. by (ii). By the definition of 371, g G P^ H^, and standard arguments show 

from this that g{C) e F(C)5© a.e. Therefore g = a.e. It follows that 971 = {0}, 
and so F is factorable by Lemma l4?2l 

Let F = G*G a.e. with G outer. We prove that ipG* € i?^0) by showing 

that ^G*H^ C Since G is outer, GH^ — iJ| for some closed subspace of (S. 

By (i), 

^PG*{GHl)=iPFHlcHl. 
Therefore il}G*H^ C H^. Suppose / £ H^^^, and consider any h £ L^. Then 

(G*/,Ml|= / {f{C),GiOh{C))^da = 0, 

because ranG(C) C 5" a.e. Thus ^G*f = a.e. It follows that ipG*H'^ C H"^^ and 
therefore VG* e i/^g). □ 

For a simple application of Theorem 14.31 suppose that is a Laurent poly- 
nomial of degree m, and choose ip to be C,™P In short order, this yields another 
proof of the operator Fejer-Riesz theorem (Theorem 12. ip . 

Another application is a theorem of Sarason [SSI P- 198], which generalizes 
the factorization of a scalar-valued function in as a product of two functions 
in (see [35l p. 56]). 

Theorem 4.4. Every G in N'^^^^ can be written G = G1G2, where Gi and G2 
belong to -/V^^igj and 

G;G2 = [G*G]^/2 and G^Gi ^ G2GI a.e. 
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Proof. Suppose first that G € H^^^y For each C G IT, write 

G(C) = f/(C)[G*(C)G(C)]^ 

where U{() is a partial isometry with initial space ran [G*(C)G(C)] 5 . It can be 
shown that U is weakly measurable. We apply Theorem 14.31 with F = [G*G]2 and 
ip = U. Conditions (i) and (ii) of Theorem 14.31 are obviously satisfied, and so we 
obtain an outer function G2 G ^£(©) such that 

G;G2 = [G*G]i/2 a.e. 
and UG^ e Set Gi = UG^ By construction Gi G if^©), 

G = U{G*G)^ = {UGl)G2 - G1G2, 

and G|Gi = G2U*UG^ = GjG^ a.e. The result follows when G e H^^^y 

The general case follows on applying what has just been shown to ip'^G, where 
is a scalar- valued outer function such that ip^G & '-' 



The standard operator generalization of Szego's theorem also follows from 
Theorem SSI 

Theorem 4.5. Let F be a weakly measurable function on T with values in £(©) 
satisfying F >0 a.e. and 

[ log+ ||i^(C)|| da < 00 and [ log+ ||F(C)"i da < 00. 
Jt Jt 

Then F is factorable. 



Proof. Since log'^ ||-F'(C)II is integrable, we can choose a scalar outer function ipi 
such that 

F^=F/\^^\' eL^(^y 

Since log"*" H^^CC)""^!! is integrable, so is log^ ll-P'i(C)^^ II- Hence there is a bounded 
scalar outer function ip such that 

We apply Theorem 14.31 to Fi with ip = ipF^^. Condition (i) is satisfied because 
ipFi — If I. Condition (ii) holds because the values of ■(/' are invertible a.e. Thus Fi 
is factorable, and hence so is F. □ 



Theorem 14.31 has a half- plane version, the scalar inner case of which is given 
in [Sni p. 117]. This has an application to the following generalization of Akhiezer's 
theorem on factoring entire functions |50[ Chapter 6] . 
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Theorem 4.6. Let F be an entire function of exponential type t, having values in 
£(25), such that F{x) > for all real x and 



Then F{x) ~ G{x)*G{x) for all real x where G is an entire function with values 
in £(©) such that exp{—iTz/2)G is of exponential type r/2 and the restriction of 
G to the upper half-plane is an outer function. 

Matrix case: 

We end this section by quoting a few results for matrix-valued functions. The 
matrix setting is more concrete, and one can do more. Statements often require 
invertibility assumptions. We give no details and leave it to the interested reader 
to consult other sources for further information. 

Our previous definitions and results transfer in an obvious way to matrix- 
valued functions. For this we choose 25 = C" for some positive integer r and 
indentify operators on C with r x r matrices. The operator norm of a matrix is 
denoted || • ||. We write L^^^,H^^ in place of ^£(0)7 -f^^^©) and || • ||oo for the 
norms on these spaces. 

Theorem 14.51 is more commonly stated in a different form for matrix- valued 
functions. 

Theorem 4.7. Suppose that F is an r xr measurable matrix-valued function having 
invertible values on T such that F > a.e. and log"*" is integrable. Then F is 
factorable if and only if log det F is integrable. 

Recall that when F is factorable, there is a unique outer G such that F = G*G 
and G(0) > 0. It makes sense to inquire about the continuity properties of the 
mapping <I> : — > G with respect to various norms. For example, see Jacob and 
Partington [37]. We cite one recent result in this area. 

Theorem 4.8 (Barclay [5]). Let F, Fn, n = 1,2, . . . , be r x r measurable matrix- 
valued functions on T having invertible values a.e. and integrable norms. Suppose 
that F — G*G and Fn — G^Gn, where G, Gn are r x r matrix-valued outer func- 
tions such that G(0) > and G„(0) > 0, n = 1, 2, . . . . Then 





if and only if 



(i) lim / ||F(C)-^^„(C)|| da^O, and 




(ii) the family of functions {logdet i^njJ^Q is uniformly integrable. 
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A family of functions {(pa}aeA ^ is uniformly integrable if for every e > 
there is a (5 > such that Jg \(pa \ da < e for all a e A whenever (j{E) < S. See [S] 
for additional references and similar results in other norms. 

A theorem of Bourgain 9 characterizes all functions on the unit circle which 
are products hg with g,h d H°°: A function f G L°° has the form f — hg where 
g^h H°° if and only if log |/| is integrable. This resolves a problem of Douglas 
and Rudin (17j . The problem is more delicate than spectral factorization; when 
1/1 = 1 a.e., the factorization cannot be achieved in general with inner functions. 
Bourgain's theorem was recently generalized to matrix-valued functions. 

Theorem 4.9 (Barclay [HE]). Suppose F e L'^^^ has invertible values a.e. Then 
F has the form F — H*G a.e. for some G,H in if and only if log \ det_F| is 

integrable. In this case, for every e > such a factorization can be found with 

||G||oo||ff|U < ll^l|oo+£. 

The proof of Theorem 14.91 in [3] is long and technical. In fact, Barclay proves 
an L^'-version of this result for all p, 1 < p < oo. 

Another type of generalization is factorization with indices. We quote one 
result to illustrate this notion. 

Theorem 4.10. Let F be an r x r matrix-valued function with rational entries. 
Assume that F has no poles on T and that det i^(C) 7^ for all ( in T. Then there 
exist integers ki < X2 < ■ ■ ■ < Xr such that 

F{z) = F_(z)diag{z-\ . . . 

where F± are r x r matrix-valued functions with rational entries such that 

(i) ^+(2) has no poles for \z\ < 1 and deti^_|_(z) ^ for \z\ < 1; 

(ii) F^{z) has no poles for \z\ > 1 including z — 00 and det_F_(z) ^ for |z| > 1 
including z — 00. 

The case in which F is nonnegative on T can be handled using the operator 
Fejer-Riesz theorem (the indices are all zero in this case). The general case is given 
in Gohberg, Goldberg, and Kaashoek [27l pp. 236-239]. This is a large subject that 
includes, for example, general theories of factorization in Bart, Gohberg, Kaashoek, 
and Ran [7J and Clancey and Gohberg [13 . 

Historical remarks: 

Historical accounts of spectral factorization appear in [51 [301 [SHI EI]- 
Briefly, the problem of factoring nonnegative matrix-valued functions on the unit 
circle rose to prominence in the prediction theory of multivariate stationary sto- 
chastic processes. The first results of this theory were announced by Zasuhin [65] 
without complete proofs; proofs were supplied by M. G. Krem in lectures. Modern 
accounts of prediction theory and matrix generalizations of Szego's theorem are 
based on fundamental papers of Helson and Lowdenslager [3T1 , and Wiener 
and Masani [63l |64] . The general case of Theorem 14.51 is due to Devinatz [15] ; 
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other proofs are given in [161 1301 150] . For an engineering view and computational 
methods, see [2H1 Chapter 8] and [55] , 

The original source for Lowdenslager's Criterion (Lemmas 12.31 and 14. 2p is 
|40j : an error in [40] was corrected by Douglas [16] . There is a generalization, 
given by Sz.-Nagy and Foias [HTJ pp. 201-203], in which the isometry may have a 
nontrivial unitary component and the shift component yields a maximal factorable 
summand. Lowdenslager's Criterion is used in the construction of canonical models 
of operators by de Branges [10] ■ See also Constantinescu _14J for an adaptation to 
Toeplitz kernels and additional references. 

5. Multivariable theory 

It is natural to wonder to what extent the results for one variable carry over to sev- 
eral variables. Various interpretations of "several variables" are possible. The most 
straightforward is to consider Laurent polynomials in complex variables zi, . . . , Zd 
that are nonnegative on the d-torus T''. The method of Schur complements in fJ3] 
suggests an approach to the factorization problem for such polynomials. Care is 
needed, however, since the first conjectures for a multivariable Fejer-Riesz theorem 
that might come to mind are false, as explained below. Multivariable generaliza- 
tions of the Fejer-Riesz theorem are thus necessarily weaker than the one- variable 
result. One difficulty has to do with degrees, and if the condition on degrees is 
relaxed, there is a neat result in the strictly positive case (Theorem 15. ip . 

By a Laurent polynomial in z — (zi, . . . , z^) we understand an expression 

mi rrid 

Q(^)= E ••• E Qk^.....kA'■■■4'■ (5.1) 

We assume that the coefficients belong to £(25), where © is a Hilbert space. With 
obvious interpretations, the scalar case is included. By an analytic polynomial with 
coefficients in £((5) we mean an analogous expression, of the form 

mi nid 

E-'-E^^-.^^^^i'^-'-^^- (5-2) 

fci=0 fcd=0 

The numbers mi, . . . ,md in ()5.ip and ()5.2p are upper bounds for the degrees of 
the polynomials in zi, . . . , Zd, which we define as the smallest values of mi, . . . , m^ 
that can be used in the representations (|5.ip and (|5.2p . 

Suppose that Q{z) has the form (|5.ip and satisfies Q{C) > for all C £ 
T*^, that is, for all ( = {C,i, ■ ■ ■ Xd) with \C^i\ = ■■■ = \C,d\ — 1- Already in the 
scalar case, one cannot always find an analytic polynomial P{z) such that (5(C) = 
P(C)*P(C), C G T'^- This was first explicitly shown by Lebow and Schreiber [55] , 
There are also difficulties in writing Q{C) = X;j=i Pj{C)*PjiC), C e T'^, for some 
finite set of analytic polynomials, at least if one requires that the degrees of the 
analytic polynomials do not exceed those of Q{z) as in the one- variable case (see 
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Naftalevich and Schreiber [Uj, Rudin [53], and Sakhnovich [5i, §3.6]). The example 
in [33] is based on a Cayley transform of a version of a real polynomial over 
called Motzkin's polynomial, which was the first explicit example of a nonnegative 
polynomial in M'', d > 1, which is not a sum of squares of polynomials. What is 
not mentioned in these sources is that if we loosen the restriction on degrees, the 
polynomial in [44] can be written as a sum of squares (see [H]). Nevertheless, for 
three or more variables, very general results of Scheiderer j57j imply that there 
exist nonnegative, but not strictly positive, polynomials which cannot be expressed 
as such finite sums regardless of degrees. 

Theorem 5.1. Let Q{z) he a Laurent polynomial in z = (zi, . . . , Zd) with coefficients 
in £(©) for some Hilbert space ©. Suppose that there is ad > such that Q{C) > dl 
for all C e T'*. Then 

r 

Q{C)^Y.P^^O*PAO, CeT'', (5.3) 

for some analytic polynomials Pi{z), . . . , Pr{z) in z = (zi, . . . , Zd) which have co- 
efficients in £(©). Furthermore, for any fixed k, the representation (I5.3P can he 
chosen such that the degree of each analytic polynomial in Zk is no more than the 
degree of Q{z) in Zk- 

The scalar case of Theorem 15.11 follows by a theorem of Schmiidgen [58] , 
which states that strictly positive polynomials over compact semialgebraic sets 
in R" (that is, sets which are expressible in terms of finitely many polynomial 
inequalities) can be written as weighted sums of squares, where the weights are 
the polynomials used to define the semialgebraic set (see also T2]); the proof is 
nonconstructive. On the other hand, the proof we sketch using Schur complements 
covers the operator- valued case, and it gives an algorithm for finding the solution. 
One can also give estimates for the degrees of the polynomials involved, though 
we have not stated these. 

We prove Theorem lS.ll for the case d — 2, following Dritschel [TS]. The general 
case is similar. The argument mimics the method of Schur complements, especially 
in its original form used in [T^. In place of Toeplitz matrices whose entries are 
operators, in the case of two variables we use Toeplitz matrices whose entries 
are themselves Toeplitz matrices. The fact that the first level Toeplitz blocks are 
infinite in size causes problems, and so we truncate these blocks to finite size. 
Then everything goes through, but instead of factoring the original polynomial 
Q{z), the result is a factorization of polynomials Q^^\z) that are close to Q{z). 
When Q(C) > SI on T'^ for some (5 > 0, there is enough wiggle room to factor Q{z) 
itself. We isolate the main steps in a lemma. 

Lemma 5.2. Let 

mi 7712 

j — — mi k— — 7rL2 
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be a Laurent polynomial with coefficients in £(25) such that Q{C) > for all 
C= (Ci,C2) m T^. Set 

j— — rni k—~'m2 

Then for each N > m2, there are analytic polynomials 

mi N 



^(-)-EE^^^ 

j=0 k=0 

with coefficients in £(25) such that 

N 



■2 I 



e = o,...,N, 



(5.4) 



(5.5) 



Proof Write 

mi ^ m2 \ mi 

Q(^)= E ( E E 

j ——m\ ^ k— — m2 j ——ni\ 

and extend all sums to run from ~oo to oo by setting Qjk — and Rj{z2) = if 
|j| > mi or I A: I > m2. Introduce a Toeplitz matrix T whose entries are the Toeplitz 
matrices Tj corresponding to the Laurent polynomials Rj{z2), that is, 

/Qjo Qj.-i Qj-2 • • -X 



/To T_i r_2 --A 



Ti To T_i 
T2 Ti To 



V 



T, 



7 



Qji Qjo Qj.-i 
Qj2 Qji Qjo 



V 



7 



j = 0, ±1, ±2, .... Notice that T is finitely banded, since Tj = for > nii. The 
identity l|2.5p has the following generalization: 

OC CXD „ 

{Th,h) = y2y2{T,_php,h,)^ / (Q(c)MC),MC))erf'T2(c). 

Here ( = (Ci,C2) and dCT2(C) = do-(Ci)c?cr(C2). Also, 



^2 ' 



MC) = EEvcfC2' 

where the coefficients are vectors in © and all but finitely many are zero, and 



h = 







^hpQ 


hi 


hp = 


hpi 


[& 




\ '■ 



P = 0,f,2, 
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It follows that T acts as a bounded operator on a suitable direct sum of copies 
of 0. Since Q(C) > on T^, T > 0. 

Fix N >m2. Set 











•■\ 






n 


n Tu 






r = 










n 










\ '■■ 




■/ 


where T- is the upper (iV + 1) x (A^ 


+- 1) block of Tj ^ 






1 


Qj-N ^ 






Qji Qjo 




Qj,-N+l 




J iV+ 1 












\QjN Qj,N- 


-1 


QjQ / 





j = 0,±l,±2,... . 

Then T' is the Toeplitz matrix corresponding to the Laurent polynomial 

mi 
j=-mi 

Moreover, T' > since it is a positive constant multiple of a compression of T. 
Thus ^(zf) > for \w\ = 1. By the operator Fejer-Riesz theorem (Theorem 12. ip . 



'^(w) = <f>(w)*$(it;), 



w 



(5.6) 



for some analytic polynomial ^{w) = X]j=?o ^j'^'' whose coefficients are {N + 1 
{N + 1) matrices with entries in £(©). Write 



where is the fc-th column in Set 



mi N 
j=Q k=0 

The identity (|5.6p is equivalent to 2mi + 1 relations for the coefficients of 'i'{w). 
The coefficients of \['(?«) are constant on diagonals, there being N + 1 — k terms 
in the fc-th diagonal above the main diagonal, and similarly below. If these terms 
are summed, the result gives 2mi -I- 1 relations equivalent to the identity 



gW(C) = i^(Cri^(C), CeT2. 



(5.7) 
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We omit the calculation, which is straightforward but laborious. To convert (|5.7p 
to the form (|5.5p . write 



Then 



F, 



(1) 



jk 
(N) 



j — 0, . . . , nil and A; = 0, . . . , A^. 



Fiz) 



Fi{z) 



\Fn{z)J 

where Fq{z), . . . , Fj^{z) are given by (|5.4p . and so (|5.7p takes the form (|5.5p . □ 



Proof of Theorem \5.1\ for the case d = 2. Suppose N > m2, and set 

™2 N + 1 

J — —'mi /c— — m2 

The values of are selfadjoint on T^, and Q{z) — Q{z) + S{z), where 

mi "12 1,1 

j — — m\ k— — 7n2 

Now choose TV large enough that ||S'(C)|| < (^,_C e Then Q(C) > on T^, and 
the result follows on applying Lemma [5?2l to Q{z). □ 

Further details can be found in [18j , and a variation on this method yielding 
good numerical results is given in Geronimo and Lai [23j . 

While, as we mentioned, there is in general little hope of finding a factoriza- 
tion of a positive trigonometric polynomial in two or more variables in terms of 
one or more analytic polynomials of the same degree, it happens that there are 
situations where the existence of such a factorization is important. In particular, 
Geronimo and Woerdeman consider this question in the context of the autoregres- 
sive filter problem [211 [25] , with the first paper addressing the scalar case and the 
second the operator-valued case, both in two variables. They show that for scalar- 
valued polynomials in this setting there exists a factorization in terms of a single 
stable (so invertible in the bidisk B^) analytic polynomial of the same degree if 
and only if a full rank condition holds for certain submatrices of the associated 
Toeplitz matrix ( [24|, Theorem 1.1.3]). The condition for operator- valued polyno- 
mials is similar, but more complicated to state. We refer the reader to the original 
papers for details. 
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Stable scalar polynomials in one variable are by definition outer, so the Geron- 
imo and Woerdeman results can be viewed as a statement about outer factoriza- 
tions in two variables. In [121 , ^ different notion of outerness is considered. As 
we saw in ij3l in one variable outer factorizations can be extracted using Schur 
complements. The same Schur complement method in two or more variables gives 
rise to a version of "outer" factorization which in general does not agree with 
that coming from stable polynomials. In [19j . this Schur complement version of 
outerness is used when considering outer factorizations for polynomials in two or 
more variables. As in the Geronimo and Woerdeman papers, it is required that 
the factorization be in terms of a single analytic polynomial of the same degree as 
the polynomial being factored. Then necessary and sufficient conditions for such 
an outer factorization under these constraints are found ( |19[ Theorem 4.1]). 

The problem of spectral factorization can also be considered in the mul- 
tivariable setting. Blower [5] has several results along these lines for bivariate 
matrix-valued functions, including a matrix analogue of Szego's theorem similar 
to Theorem 14.71 His results are based on a two- variable matrix version of Theo- 
rem 15.11 and the arguments he gives coupled with Theorem 15.11 can be used to 
extend these results to polynomials in d > 2 variables as well. 

6. Noncommutative factorization 

We now present some noncommutative interpretations of the notion of "several 
variables," starting with the one most frequently considered, and for which there 
is an analogue of the Fejer-Riesz theorem. It is due to Scott McCuUough and comes 
very close to the one- variable result. Further generalizations have been obtained 
by Helton, McCuUough and Putinar in [33]. For a broad overview of the area, 
two nice survey articles have recently appeared by Helton and Putinar [34] and 
Schmiidgen [59] covering noncommutative real algebraic geometry, of which the 
noncommutative analogues of the Fejer-Riesz theorem are one aspect. 

In keeping with the assumptions made in [l^ , all Hilbert spaces in this section 
are taken to be separable. Fix Hilbert spaces 6 and Sj, and assume that is infinite 
dimensional. 

Let S be the free semigroup with generators ai, . . . , a^. Thus S is the set of 
words 

w ^ aj^ ■ ■ -aj^, G {l,...,(i}, fc = 0,l,2,..., (6.1) 

with the binary operation concatenation. The empty word is denoted e. The length 
of the word (|6.ip is \w\ = k (so |e| = 0). Let Sm be the set of all words (|6.ip of 
length at most m. The cardinality of Sm is im — ^ + d + (P + ■ ■ ■ + (f^ . 

We extend 5 to a free group G. We can think of the elements of G as words 
in oi, . . . , ad, , ■ ■ ■ , o^^, with two such words identified if one can be obtained 
from the other by cancelling adjacent terms of the form aj and aj^. The binary 
operation in G is also concatenation. Words in G of the form h =^ v^^w with 
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v,w S play a special role and are called hereditary. Notice that a hereditary 
word h has many representations h — v^-^w with w,w e 5. Let Hm be the set of 
hereditary words h which have at least one representation in the form h — v^^w 
with v,w € Sm- 

We can now introduce the noncommutative analogues of Laurent and analytic 
polynomials. A hereditary polynomial is a formal expression 

g= ^ h®Qh, (6.2) 

where Qu G for all h. Analytic polynomials are hereditary polynomials of 

the special form 

P= w(E)P^, (6.3) 

where Pw G S,{<3) for all w. The identity 

Q = P*P 

is defined to mean that 

Next we give meaning to the expressions Q{U) and P{U) for hereditary and 
analytic polynomials (|6.2[) and (|6.3[) and any tuple U = {Ui, . . . , Ud) of unitary 
operators on i^. First define G for any w G S hy writing w in the form 

(|6.ip and setting 

By convention, C7^ = / is the identity operator on i5. If ft, G © is a hereditary word, 
set 

jjh ^ {jjvyijw 

for any representation h = v~^w with w, w G S*; this definition does not depend on 
the choice of representation. Finally, define Q{U), P{U) G ® ©) by 

Q(C/)= U''®Qh, P{U)^ Yl U'^^P^. 

The reader is referred to, for example, Murphy [43l §6.3] for the construction of 
tensor products of Hilbert spaces and algebras, or Palmer, [331 §1.10] for a more 
detailed account. 
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Theorem 6.1 (McCuUough 42J). Let 

Q = ^ h®Qh 

he a hereditary polynomial with coefficients in £(©) such that Q(U)) > for every 
tuple U = (J7i, . . . , Ud) of unitary operators on ^j. Then for some I < f-m, there 
exist analytic polynomials 

Pj = ^ w® Pj^^, j = !,...,£, 
with coefficients in £(6) such that 

Q = p*Pi + --- + p;Pi. 

Moreover, for any tuple U = {Ui, . . . , Ud) of unitary operators on ^j, 

Q{u) = Pi{uyPi{u) + ■■■ + Piiuypeiu). 

In these statements, when © is infinite dimensional, we can choose £ — I. 

As noted by McCuUough, when d = 1, Theorem 16.11 gives a weaker version 
of Theorem 2.1. However, Theorem 2.1 can be deduced from this by a judicious 
use of Beurhng's theorem and an inner-outer factorization. 

McCuUough's theorem uses one of many possible choices of noncommutative 
spaces on which some form of trigonometric polynomials can be defined. We place 
this, along with the commutative versions, within a general framework, which we 
now explain. 

The complex scalar-valued trigonometric polynomials in d variables form a 
unital ^-algebra ^P, the involution taking z" to z~^\ where for n = (ni, . . . ,nd), 
—n = (— ni, . . . , — nrf). If instead the coefficients are in the algebra £ — £(©) for 
some Hilbert space ©, then the unital involutive algebra of trigonometric polyno- 
mials with coefficients in £ is £. The unit is 101. A representation of^ ®)€ 
is a unital algebra *-homomorphism from *p 05 £ into for a Hilbert space Sj. 

The key thing here is that zi, . . . , Zd generate £, and so assuming we do not mess 
with the coefficient space, a representation tt is determined by specifying 7r(zfc), 
k = I, . . . ,d. 

First note that since z^z^ = 1, 7r(zfc) is isometric, and since = z'^^, we 
then have that Tr{zk) is unitary. Assuming the variables commute, the ZkS gener- 
ate a commutative group G which we can identify with Z'^ under addition, and 
the irreducible representations of commutative groups are one dimensional. This 
essentially follows from the spectral theory for normal operators (see, for exam- 
ple, Edwards [21] p. 718]). However, the one-dimensional representations are point 
evaluations on T*^. Discrete groups with the discrete topology are examples of 
locally compact groups. Group representations of locally compact groups extend 
naturally to the algebraic group algebra, which in this case is and then on to 
the algebra (8) £ by tensoring with the identity representation of £. So a seem- 
ingly more complex way of stating that a commutative trigonometric polynomial 
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P in several variables is positive / strictly positive is to say that for each (topo- 
logically) irreducible unitary representation tt of G, the extension of tt to a unital 
*-representation of the algebra ^P^C, also called tt, has the property that 7r(P) > 
/ 7r(P) > 0. By the way, since T'' is compact, 7r(P) > implies the existence of 
some e > such that n{P - el ® 1) = 7r(P) - el > 0. 

What is gained through this perspective is that we may now define noncom- 
mutative trigonometric polynomials over a finitely generated discrete (so locally 
compact) group G in precisely the same manner. These are the elements of the 
algebraic group algebra ^ generated by G; that is, formal complex linear com- 
binations of elements of G endowed with pointwise addition and a convolution 
product (see Palmer [45l section 1.9]). Then a trigonometric polynomial in *p ® £ 
is formally a finite sum over G of the form P — J2g 9 ® Pg where Pg E C for all g. 

We also introduce an involution by setting g* = g~^ for g G G. A trigonomet- 
ric polynomial P is self adjoint if for all g, Pg* = P* . There is an order structure 
on selfadjoint elements defined by saying that a selfadjoint polynomial P is posi- 
tive / strictly positive if for every irreducible unital *-representation tt of G, the 
extension as above of tt to the algebra £ (again called tt), satisfies 7r(P) > / 
7r(P) > 0; where by 7r(P) > we mean that there exists some e > independent 
of TT such that 7r(P — e(l ® 1)) > 0. Letting fl represent the set of such irreducible 
representations, we can in a manner suggestive of the Gel'fand transform define 
P{tt) = 7r(P), and in this way think of as a sort of noncommutative space on 
which our polynomial is defined. The Gel'fand-Raikov theorem (see, for example. 
Palmer [551 Theorem 12.4.6]) ensures the existence of sufficiently many irreducible 
representations to separate G, so in particular, ^ 0. 

For a finitely generated discrete group G with generators {oi, . . . ,ad}, let 
5 be a fixed unital subsemigroup of G containing the generators. The most in- 
teresting case is when S is the subsemigroup generated by e (the group identity) 
and {fli, . . . , a^}- As an example of this, if G is the noncommutative free group 
in d generators, then the unital subsemigroup generated by {ai, . . . ,ad} consists 
of group elements w of the form e (for the empty word) and those which are an 
arbitrary finite product of positive powers of the generators, as in (|6.ip . 

We also need to address the issue of what should play the role of Lau- 
rent and analytic trigonometric polynomials in the noncommutative setting. The 
hereditary trigonometric polynomials are defined as those polynomials of the form 
P = J2j '^ji^j^ Pjj where Wji,Wj2 G S. We think of these as the Laurent poly- 
nomials. Trigonometric polynomials over S are referred to as analytic polynomials. 
The square of an analytic polynomial Q is the hereditary trigonometric polynomial 
Q*Q. Squares are easily seen to be positive. As a weak analogue of the Fejer-Riesz 
theorem, we prove a partial converse below. 

We refer to those hereditary polynomials which are selfadjoint as real hered- 
itary polynomials, and denote the set of such polynomials by H. While these 
polynomials do not form an algebra, they are clearly a vector space. Those which 
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are finite sums of squares form a cone C in H (that is, C is closed under sums 
and positive scalar multiplication) . Any real polynomial is the sum of terms of the 
form 1 (8) A or W2'Wi B + wlw2 ^ B* , where wi,W2 ^ S and A is selfadjoint. The 
first of these is obviously the difference of squares. Using w*w = 1 for any w e G, 
we also have 

^2^1 <»B + WiW2 (g) B* = (Wl ® B + W2 ® l)*(wi ® B + ^2 «) 1) - 1 ® (1 + B*B). 
Hence H ^C-C. 

For A,B <E £(^) and wi,W2 e S, 

Q<{wi®A + w2® B)*{wi ®A + W2®B) 
< {wi(®A + w2(E)B)*iwi(g)A + w2(E)B) 

+ {wi (g> A ~ W2 (E) B)*{wi (g) A~W2'S> B) 
= 2{1®A*A + 1(SB*B) 
<(||A|p + ||Bf)(l0l). 

Applying this iteratively, we see that for any P G H, there is some constant 
< a < oo such that al ± P G C. In other words, the cone C is airchimedean. 
In particular, 1 ® 1 is in the algebraic interior of C, meaning that if P G _ff , then 
there is some < <o < 1 such that for all < < < to, t{l ® 1) + (1 - t)P e C. 

Theorem 6.2. Let G be a finitely generated discrete group, P a strictly positive 
trigonometric polynomial over G with coefficients in Then P is a sum of 

squares of analytic polynomials. 

Proof. The proof uses a standard GNS construction and separation argument. 
Suppose that for some e > 0, P — e(l 1) > but that P ^ C. Since C has 
nonempty algebraic interior, it follows from the Edelheit-Kakutani theoren{3 that 
there is a nonconstant linear functional X : H ^ such that A(C) > and 
A(P) < 0. Since A is nonzero, there is some R G H with A(P) > 0, and so since 
the cone C is archimedean, there exists a > such that a{l (g) 1) — P G C. From 
this we see that A(l (E) 1) > 0, and so by scaling, we may assume A(l ® 1) = 1. 
We next define a nontrivial scalar product on H by setting 
{wi (g)A,W2(E>B) = A(u;2 wi ® B*A) 

and extending linearly to all of H. It is easily checked that this satisfies all of the 
properties of an inner product, except that (w ® v4, w A) =0 may not necessar- 
ily imply that w^A = 0. Even so, such scalar products satisfy the Cauchy-Schwarz 
inequality, and so N ~ {w (S> A : {w A, w ® A) = 0} is a vector subspace of H. 
Therefore this scalar product induces an inner product on H/N, and the comple- 
tion of H/N with respect to the associated norm makes H/N into a Hilbert 
space. 

^ 1361 Holmes, Corollary, §4B. Let A and B be nonempty convex subsets of X, and assume 
the algebraic interior of A, cor(j4) is nonempty. Tlien A and B can be separated if and only if 
cot{A) n B = 0. 
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We next define a representation tt : H £,{9j) by the left regular represen- 
tation; that is, Tr{P)[w A] = [P{w ® A)], where [•] indicates an equivalence class 
in H/N. Since P > e(l (g) 1) > for some e > 0, P - e/2(l ® 1) > 0. Suppose that 
P^C. Then 

A((P - e/2(l ® 1)) + e/2(l (g, 1)) = A(P - e/2(l (g) 1)) + e/2 < 0. 



and so 7r(P — e/2(l ® 1)) ^ 0. The representation tt obviously induces a uni- 
tary representation of G via n(ai) = 7r(ai g) 1), where Ui is a generator of G. 
The (irreducible) representations of G are in bijective correspondence with the 
essential unital ^-representations of the group C*-algebra C*{G) (Palmer, [151 
Theorem 12.4.1]), which then restrict back to representations of H. Since unitary 
representations of G are direct integrals of irreducible unitary representations (see, 
for example, Palmer, [JSl p. 1386]), there is an irreducible unitary representation 
tt' of G such that the corresponding representation of H has the property that 
7r'(P — e/2(l ® 1)) 2 0, giving a contradiction. □ 

The above could equally well have been derived using any C*-algebra in the 
place of One could also further generalize to non-discrete locally compact 

groups, replacing the trigonometric polynomials by functions of compact support. 

We obtain Theorem 15.11 as a corollary if we take G to be the free group in 
d commuting letters. On the other hand, if G is the noncommutative free group 
on d letters, it is again straightforward to specify the irreducible representations 
of G. These take the generators (ai, . . . , ad) to irreducible d-tuples {Ui, . . . , Ud) of 
(noncommuting) unitary operators, yielding a weak form of McCuUough's theorem. 

As mentioned earlier, it is known by results of Scheiderer [57] that when G is 
the free group in d commuting letters, c? > 3, there are positive polynomials which 
cannot be expressed as sums of squares of analytic polynomials, so no statement 
along the lines of Theorem 16.11 can be true for trigonometric polynomials if it is 
to hold for all finitely generated discrete groups. Just what can be said in various 
special cases is still largely unexplored. 



Appendix A. Schur complements 

We prove the existence and uniqueness of Schur complements for Hilbert space 
operators as required in Definition 13.11 

Lemma A.l. Let T e where Sj is a Hilbert space. Let A be a closed subspace 

of S), and write 



Hence 



(7r(P - e/2(l ® 1))[1 ® 1], [1 ® 1]) < -e/2. 
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Then T > if and only if A > 0, C > 0, and B — C'^GA'^ for some contraction 
G G ) . The operator G can be chosen so that it maps ran A into ran C and 

is zero on the orthogonal complement of ran A, and then it is unique. 

Proof. U B = C^GA^ where G G £(^, is a contraction, then 



,C3G C5(/-GG*)3y V (/-GG*)3C5 
Conversely, if T > 0, it is trivial that A>0 and G > 0. Set 

Then ^ = NiN^ and G = and so there exist partial isometries Vi G ij) 

and V2 G with initial spaces ran^ and ranG such that = V1A2 and 

TV* = V2Ci. Thus B = iVsiVi* = G^GA^, where G = ^3*^1 is a contraction. 
By construction G has the properties in the last statement, and clearly such an 
operator is unique. □ 

Lemma A. 2. Let be a Hilbert space, and suppose T G ^{^), T > 0. Let ^ be a 
closed subspace of S), and write 

^"{b ^ 

Then there is a largest operator S > in £,{A) such that 

'A~S B 



B C 



> 0. (A.l) 



It is given by S = A2{I- G*G)A2, where G G £(J?, j?^) is a contraction which 
maps ran A into ranG and is zero on the orthogonal complement of ran A. 

Proof. By Lemma lATl we may define S — A^{I — G*G)A^ with G as in the last 
statement of the lemma. Then 



A-S B*\ _fA2G*GA2 A2G*C2\ _ fA2G' 



{ga^ G^) > 0. 



Consider any X > in £(^) such that 

'A-X B* 
B C 



> 0. 



Since A > X > 0, we can write X = A2KA2 where K G and < X < /. 

We can choose K so that it maps ran A into itself and is zero on (ran A)-'-. Then 

A-X B*\_fA-AiKAi A^G*Ci 
B C ) - [ C^GA^ G 

A^ \ fl -K G*\ (A^ 
G5 i V G / j I C"2 



30 



Michael A. Dritschel and James Rovnyak 



By our choices G and K, we deduce that 

By Lemma G = Gi{I — K)^ where Gi G £(^, is a contraction. Therefore 
G*G<I-K, and so 

X = yl^i^A^ <A^I- G*G)A^ = S. 

This shows S is maximal with respect to the property (jA.ip . □ 
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